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Abstract. It is generally accepted that the PSR B 1509- 
58 is associated with the supernova remnant (SNR) 
MSH 15-52 (G 320.4-01.2). The spin-down age of the pul- 
sar is ~ 1700 years, while the size and the general appear- 
ance of the SNR suggest that this system is much older. 
A few possible explanations of this discrepancy have been 
put forward. We offer an alternative one and suggest that 
the high spin-down rate of the pulsar characterizes only 
a relatively short period of its (present) spin history, and 
that the enhanced braking torque is connected with the 
interaction between the pulsar's magnetosphere and the 
dense matter of a circumstellar clump (created during the 
. late evolutionary stages of the supernova (SN) progeni- 
tor star). Our suggestion implies that the "true" age of 
PSR B 1509-58 could be much larger than the spin-down 
age, and therefore the SNR MSH 15-52 is a middle-aged 
(y-j , remnant similar to the Vela SNR (G 263.9-3.3). We also 
suggest that the dense (neutral) gas of the circumstellar 
clump could be responsible for the enhanced neutral hy- 
■ drogen absorption towards PSR B 1509-58, and that the 
optical emission of an optical counterpart for PSR B 1509- 
58 should rather be attributed to a bow shock around this 
i; pulsar than to the pulsar itself. 
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1. Introduction 



PSR 



B 1509-58 
Manchester et al. 1982D 



dSeward fc Harnden 1982 , 
is situated not far from the 
geometrical centre of the SNR MSH 15-52 (or G 320.4- 
01.2), and their association is beyond any doubt. However, 
this association causes a number of difficulties for the 
understanding and interpretation of the observational 
data. The problem is that the size and the general 
appearance of the SNR suggest that it should be much 
older than it follows from the pulsar age estimates. 

It is usually assumed that the rotational frequency 
VL of a. pulsar decreases according to the relation Cl = 
— K^}", where K depends upon the physics of the slow- 
down mechanism, and n ~ flfl/fl^ is the braking in- 
dex. Assuming constant K and n, and provided that the 
initial spin period P; of the pulsar was much smaller 
than the current period P — 27r/ri, one can estimate 
the characteristic spin-down age r — P/{n — l)P. For 
P ~ 0.15 s, P ~ 1.49 • 10-12 ss^i and n = 2.84 



(Weisskopf et al. 1983, Manchester et al. 1985, Kaspi et 
al. |1994| ), one derives an age of PSR B 1509-58 of ~ 
1700 years, i.e. it is nearly as young as the Crab pul- 
sar. The spin-down age could be even less by a factor 
of [1 - {Pi/P)"^'^] -\ if the pulsar was b orn with P; 
of ~ 0.1 s (see e.g. 3pruit fc Phinney 1998 ). These esti- 
mates are at odds with t he age estimates for MSH 15- 
52 (Seward et al. 1983, van den Berg fc Kamper 1984 , 
Kamper et al. 1995 ), which show that the SNR is a much 
older object. 

To reconcile the ages of the pulsar and the SNR, Se- 
ward et al. ( 1983| ) considered two possibilities: 1) MSH 15- 
52 is a young SNR, and 2) PSR B 1509-58 is an old 
pulsar. The first one implies (in the framework of the 
Sedov- Taylor model) that the SN explosion was very 
energetic and occu red in a tenuous medium (see also 
Bhattacharya 1990| ). This point of view is generally ac- 
cepted (e.g. Kaspi e t al. 1994 , Greiveldinger et al. 1995 , 
Trussoni et al. 19"96| , Gaenzlcr et al. 1999). The second 
possibility implies that r is at least few times shorter than 



the "true" age (Seward et al. 1983|). This poss ibility was 
re-examinated by p31andford fc Romani (1988) . Assuming 
that the pulsar spin-down is mostly due to the electromag- 
netic torque, they suggested that the torque grew within 
the last ~ 10'^ years due to the growth of the pulsar's mag- 



netic field (see also Muslimov & Page 1996). In this case, 



the coefficient K is an increasing function of time and 
therefore the "true" age of the pulsar could be as large as 
follows from the age estimates for the SNR. In this paper 
we offer an alternative explanation for the increase of the 
braking torque (Sect. 2), viz. we suggest that it could be 
episodically enhanced due to the interaction of the pul- 
sar's magnetosphere with dense clumps of circumstellar 
matter (Sect. 3). Sect. 4 deals with some issues related to 
our suggestion. 



2. The spin-down of PSR B 1509-58 

It is known that the electromagnetic torque acting on 
a rotating, magnetized body (e.g. a neutron star) im- 
mersed in a plasma is enhanced as compared with the 

see also Istomin 



torque in vacuum (Ginzburg 1971 



1994) 



It was mentioned by Istomin (1994) that for the in- 
crease of the slow-down torque of a pulsar it is suffient 
to have a dense plasma in the vicinity of the light cylin- 
der, since just in this region the pulsar loses its rotational 
energy due to the acceleration of particles. The parti- 
cles of the ambient medium penetrating into the region 
of the light cylinder are accelerated there to velocities 
comparable with the speed of light and then leave this 
region (Istomin 1994). This presumably equatorial out- 
flow (cf. Brinkmann et al. 1985, King fc Cominsky 1994) 



carries away the pulsar's angular momentum and is re- 
sponsible for the enhanced braking of the pulsar. We sug- 
gest (see also Gvaramadze 1999a; cf. Yusifov et al. 1995 



[stomin fc Romberg 2000) that just this effect is responsi- 
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ble for the present high spin-down rate of PSR B 1509-58, 
i.e. that the pulsar loses its rotational energy mainly due 
to the acceleration of protons of the ambient medium ar- 
riving at the light surface at the rate M: 

\E\ = \mn\ = 7pMc2 , (1) 
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where / ~ lO^^gcm^ is the moment of inertia of the 
pulsar, c is the speed of light, and 7p is the Lorentz 
factor of accelerated protons. For H. — 41.89 s^^ and 
n = -4. 25-10 '^" s^2 (g_g^ K:j^gpj |1994| ), and 7p ~ 1 
(Istomin |l99^, one has M ~ 2.2 • 10^'' gs~^ We also sug- 
gest that the pulsar moves through the inhomogeneous 
ambient medium and episodically plunges into clumps of 
dense matter. If the density of clumps is sufficiently high 
(see Sect. 3), one can expect a temporary increase of f2 
every time the pulsar travels through a clump (we believe 
that just this situation takes place now). The correspond- 
ing spin-down age will be less than estimated when the 
pulsar moves through the low-density interclump medium. 

It is clear that the presence of radio emission of the 
pulsar means that the ambient medium does not pene- 
trate far beyond the light surface. Assuming that the ram 
pressure of the accreting mediu m is equal to the magneti c 
pressure at the light surface (cf . [King fc Cominsky 1994 ) , 



one has an estimate of the surface magnetic field of the 
pulsar^ : 

i?* 



7/2 , 



1/2 



2.9-10^^0 



L /rt 

where G is the gravitational constant. A/* — 1.4 Mq and 
= 10^ cm are the mass and the radius of the pulsar, and 
tl — c/il ~ 7.1 • 10^ cm is the radius of the light cylinder; 
for simplicity we assumed that the pulsar magnetic field 
is dipolar. Given this value of _B», one can estimate the 
"vacuum" values of f2 and r: = -2ii'^n^ /3c^I ~ -1.36- 
10"^"'^ s~^, where /i — B^,r^ is the magnetic momentum of 
the pulsar, and To < -Vl/2tlo ~ 4.9 • lO'^yr ~ 30t (here 
we assumed that the "vacuum" braking index is equal to 
3). This means that the "true" age of the pulsar could be 
as large as tq (provided that Pi ^ P) and that the SNR 
MSH 15-52 could be a middle-aged remnant similar to the 
Vela SNR (G 263.9-3.3). 



3. Dense clumps in the central part of SNR 
MSH 15-52 

We suggest that PSR B 1509-58 and SNR MSH 15-52 
are th e remnants of the SN explosion of a massive 
star ( iGvaramadze 1999b ). In this case, the struc- 
ture of MSH 15-52 could be determined by the 
interaction of the SN blast wave with the ambi- 
ent medium reprocessed by the joint action of the 
ionizing emission and stellar wind of the SN pro- 



gcnitor star ( 


McKee et al. 1984, 


ShuU ct al. 1985, 


Ciotti & D'Ercole 1989, 


Chevalier & Liang 1989, 


Franco ct al. 1991 


D'Ercole 1992, Gvaramadze 1999b 


c. 


2000a 


). The outer shell of the SNR could arise due 


to 



the abrupt deceleration of the SN blast wave after it 
encounters the density jump at the edge of the bubble 
created by the fast stellar wind during the main-sequence 



^ Note that this estimate is about ten times smaller 
t han wliat follows from the usual dipole formula 
([Weisskopf et al. 1983|). 



or the Wolf-Rayet (WR) stages. On the other hand, 
some structures in the central part of the SNR|| could 
be attributed to the interaction of the SN blast wave 
with the circumstellar material lost during the late 
evolutionary stages of the SN progenitor star (this is 
the material that determines the appearance of young 
type I I SNRs, e.g. SN 1987A (e. g. jMcCray 1993D or Gas A 
(e.g. Garcia-Segura et al. 1996| , Borkowski et al. 1996D ). 
During the red supergiant (RSG) stage a massive star 
loses a significant part (about two thirds) of its mass in 
the form of a slow, dense wind. This matter occupies 
a compact region with a characteristic radius of few 
parsecs (the high-pressure gas in the main-sequence 
bubble significantly affects the spr ea ding of this r e- 
gion ( IChcvalicr fc Emmering 1989i [D'Ercole 199^ ). 
Before the SN exploded, the progenitor star (of mass 
> 15 — 20 Mq) becomes fo r a short time a WR star 
(e.g. Vanbeveren et al. 1998| ). At this stage, the fast 



stellar wind sweeps up the slow RSG wind and creates 
a low-density cavity surrounded by a shell of swept-up 
circumstellar matter. The shell expands with a nearly 
constant velocity v^h ^ (Mwr^^wr"RSg/3Mrsg)^^^, 
where Mwr, Mrsq and wwR, 'i^RSG are, correspondingly, 
the mass-loss rates and wind velocities during the WR 



and RSG stages (e.g. Dyson 1981), until it catches up the 
shell separating the RSG wind from the main-sequence 
bubble. For parameters typical for RSG and WR winds, 
one has Wgh — 100 — 2G0kms~^. The interaction of two 
circumstellar shells results in Rayleigh- Taylor and other 
dynamical instabilities, whose development is accom- 
panied by the formation of dense clumps moving wit h 
radial velocities of ~ Vsh ( Garcia-Segura et al. 1996 ). 
The dense clumps could originate much closer to the 
SN progenitor star due to the stellar wind acceleration 
during the transition from the RSG to the WR stage 
( Brighenti fc D'Ercole 1997 ). The number density of 
clumps is > 10^ cm~'^ provided they are not fully ionized 
and were able to cool to a temperature of < 10^ K 
( Brighenti fc D'Ercole 1997 ). Direct evidence of the 
existence of high-density clumps close to the SN explosion 
sites follows from observations of young SNRs. For exam- 
ple, the optically emitting gas of quasi-stationary flocculi 

10^ cm-3 



in Gas A is characterized by a density of 
and a temperature of > IC* K (e.g. Lozinskaya 1992 ). 
Assuming that the optical emission of a floccule comes 
from an ionized "atmosphere" around the neutral core, 
one can estimate the density of the core to be > 10^ cm~^, 
provided that the temperature of the core is < 10^ K. 
Similar estimates could also be derived from observations 
of the optical ring around SN 1987A, the inner ionized 
"skin" of which has nearly the same parameters (e.g. 



^ Note that these structures could be significantly offset 
from the geometrical centre of a middle-aged SNR due to the 



proper motion of the SN progenitor star (Gvaramadze 1999b 
2000a| Jb|). 
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Plait et al. 1995) as the optically emitting gas of flocculi 
in CasA, or from observations of some other young 



the spin-down luminosity \E\ is transferred to the am- 
bient medium (cf. Kochanek 1993). This condition can 



SNRs (e.g. phugai 1993^ [Chugai fc Danziger 1994| ). The 
radial velocity of flocculi in CasA ranges from ~ 80 to 
~ 400kms~^ (e.g. |Lozinskaya 19921 ). 

Initially, the new-born pulsar moves through the low- 
density cavity created by the fast wind of the presuper- 
nova star until it plunges into the first dense clump on 



be re-written as a < {'-fp/4)^'-v^/c ~ 10^ 7" w*,5o (cf. 
Kochanek 1993| , [Manchester et al. 1995 ), i.e. a should be 
much smaller than the usually adopted value of ~ 1 (e.g. 
Kulkarni & Hester 19881 ICordes et al. 1993| ). Weak cou- 



pling {a <^ 1) of the pulsar wind with the ambient medium 
is consistent with an outflow composed of highly relativis- 



its way. This happens at the moment t 



where tic particles (e.g. Kochanek 1993 and references therein) 



1 — 2 pc is the radius of the cavity, — Vp 



Vp and Vci are respectively the velocities of the pulsar and 
the clump. For Vp ~ 150kms~^ and Wci — 100 kms^^, one 
has t ~ (2 — 4) • 10* years^. Let us assume that all mat- 
ter captured inside the accretion radius race = 2GM^/vl 
(w* 3> Cs, where Cg is the sound speed in the cold, dense cir- 
cumstellar clump) of the pulsar moving through the clump 
penetrates into the region of the light cylinder, where it is 
accelerated to relativistic velocities and then leaves this re- 



gion in the form of equatorial outflow (cf. Istomin 1994 and 
King & Cominsky 1994). The rate at which the ambient 
medium arrives at the light cylinder could be estimated 

as 



M 



cirripV^ 



(2) 



where rid is the number density of the clump, and nip is 
the mass of a proton. If the pulsar braking is indeed mainly 
due to the acceleration of circumstellar protons arriving at 
the light surface, then one obtains from (P and (||) that 
rici ~ 3.8 • 10^ 50 cm~^, where w*.5o = u»/50kms~^. 
The time it takes for PSR B 1509-58 to cross the clump, 
t ~ Ici/v*, where Id is the characteristic size of clumps, 
should be larger than the time since the pulsar discovery, 
i.e. t > 30 years. This requirement results in la > 5 ■ 
10i5w*,50 cm and M^i > 0.0002ti3 ggM©, where M^i is the 
characteristic mass of clumps. For the SN progenitor star 
of mass > 15Mq, the mass of the circumstellar gas (i.e. 
the matter lost during the RSG stage) is about IOMq, and 
the number of clumps is < 5 • 10*. For the current radius 
of the region occupied by circumstellar matter of about 
10 pc ( iGvaramadze 1999b ), one finds a covering factor of 
the clumpy circumstellar material (i.e. the fraction of the 
sphere occupied by clumps) ~ 10~*. 

For accretion to occur, the standoff radius rg of the 
bow shock (formed by the outflow of relativistic particles) 
should be less than race- For the spherically symmetric 
outflow, one has rg = (a|i?|/47rnciTOpCD^)*/^ < race = 
2GA'h/vl, where we assume that only a fraction a < 1 of 



No prop er mo tion for PSR B 1509- 58 has yet been detected 
(Kaspi et al. L994, Gaensler et al. 1999), though the recent find- 
ing of the c andidat e optical counterpart of this pulsar (Caraveo 
et al. L994, see also Mignami et al. 1998 and Sect. 4) allows one 



to hope that it will be detected in due course. For t he distance 
to the pulsar of 5^5 kpc (e.g. Gaensler et al. 1999| ), a velocity 
of 150kms~^ corresponds to a proper motion of ~ 6masyr~^; 
this is much smaller than the upper limit s on the pulsar's 
pro per nr otion given by Kaspi et al. (1994) and Gaensler et 
al. ([L999|). 



Alternatively, if the outflow of relativistic particles is con- 
fined to the vicinity of the rotational equatorial plane, one 
can expect that the ambient matter accretes onto the pul- 
sar's magnetosphere along the polar directions. Another 
possibility is that the ambient matter penetrates in the 
pulsar wind bubble through instabilities in the bow shock 
front. In the latter both cases rg could be larger than race, 
and one can adopt a ~ 1 (see next section). 

4. Discussion 

We now discuss some consequences of our proposal that 
the braking of PSR B 1509-58 is mostly due to the interac- 
tion of the pulsar's magnetosphere with the dense matter 
of a circumstellar clump. 

First, we consider the contribution of the circumstel- 
lar matter to the neutral hydrogen absorption toward the 
pulsar. The low covering factor of clumps (see Sect. 3) 
implies that this contribution is 



ns{r)dr + nci(?'ci - ^-acc) 



(3) 



where ng (r) = M / ATrnipr'^ vs is the number density of the 
infalling gas captured inside the accretion radius, vs = 
(2G'M*/r)^/^ is the free-fall velocity, is the line of sight 
thickness of the dense neutral gas located between the 
pulsar and the observer (rd < Zci)- 



It was pointed out by Strom (|l994| ) that a ROSAT 
observation of the SNR MSH 15-52 indicates a larger 
absorption towards PSR B 1509- 58 than seen from 
the bright northwest part of the SNR's shell (known as 
RCW89). The comparison of neutral hydrogen absorption 
data (see Grcivcldingcr et al. 1995, [Trussoni et al. 1996 



Tamura et al. 1996 



Marsden et al. 1997 



Rots et al. 1998) shows that the excess of absorption to- 



wards the pulsar could be as large as ~ (1 — 5) • 10 cm~ . 
This discrepancy could be interpreted as a sign that the 
pulsar is more distant than the SNR, and therefore that 
these two objects are not physically associated with each 



other (Strom 1994). It is also possible that "the spectral 
analysis is not detailed enough to provide the correct 
parameters" ( Trussoni et al. 1996[ ). Another possibility 
is that the HI column density distribution is really 
inhomogeneous across the SNR (cf. [Trussoni et al. 1996 ). 
We favour the last possibility and suggest that the excess 
of absorption towards PSR B 1509-58 is due to the dense 
neutral gas around the pulsar. One can use Eq. (^ to set 
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an upper limit on rd. For the parameters adopted above, 
and assuming that A^ci < (1 ~ 5) 
rci < (0.2-1.2).10iX5o 



cm. This estimate shows that 
if our explanation of the age discrepancy is correct, then 
one might expect that in the near future (i.e. after a lapse 
of ~ J'ci/f*,!!, where || is the line of sight component 
of u*) the first derivative of the pulsar's spin period will 
suffer a significant decrease. 

Second, let us discuss the candidate optical counter- 



part for PSR B 1509-58 proposed by Caraveo et al. (1994) 



Caraveo et al. pointed out that the luminosity of the op- 
tical counterpart (Ly ^ 6.5 • lO^^dg ergss"^) exceeds by 
a few orders of magnitude the value derived for magne - 
tospheric optical emission of young pulsars (Pacini 1971 ) . 
This fact together with the negative result of searching 
of optical pulsations at the radio period led to the con- 
clusion that the proposed identification could be erro- 
neous (Mignami et al. 199S ; see also Shearer et al. 1998 



Chakrabarti fc Kaspi 1998 ). We suggest, however, that 



the observed optical emission should rather be attributed 
to the bow shock around the pulsar than to the pulsar 
itself. This suggestion is supported by the estimate of the 
total luminosity of the bow shock, Lt — S'ncif*(mpwJ/2), 
where S ~ 7rr^ and rb ~ S^/^rg are the area and the char- 
acteristic radius of the bow shock, respectively. For the 
adopted parameters^ this gives Lt — a ■ 10'^'^ ergs s~^, i.e. 
Ly ~ Lt for a ~ 1. It is obvious that if our suggestion is 
correct, there can be no correspondence between the ob- 
served luminosity and the luminosity expected from the 



results of Pacini (1971). The optical pulsations should be 
absent as well. 

To conclude, we point out a curious coincidence 
of the accretion rate derived in Sect. 2 with accre- 
tion rates required in accretion-based models to explain 
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